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Abstract 

This is the first in a series of articies on the numerical solution of 
Friedrich's conformal field equations for Einstein's theory of gravity. We 
will discuss in this paper why one should be interested in applying the 
conformal method to physical problems and why there is good hope that 
this might even be a good idea from the numerical point of view. We 
describe in detail the derivation of the conformal field equations in the 
spinor formalism which we use for the implementation of the equations, 
and present all the equations as a reference for future work. Finally, we 
discuss the implications of the assumptions of a continuous symmetry. 

1 Introduction 

Much of the current work in numerical and experimental relativity is devoted 
to obtaining information about the gravitational radiation that is emitted by 
astrophysical processes which are taking place in our universe. The goal is to 
obtain wave forms, i.e., the "finger prints" by which different processes can be 
identified when the gravitational waves are detected by laser interferometers like 
LIGO or VIRGO. The common way to describe such a system within Einstein's 
theory of gravity is by way of an idealization where the system is considered 
as being so far away from the rest of the universe that the influence of the 
latter can be neglected (cf. (l) for a clear discussion of what is involved in this 
process). Then, intuitively, the fields far away from the source should decay 
so that the space-time becomes asymptotically flat. The detectors are then 
idealized as observers which are located "at infinity" where they can gather the 
gravitational radiation coming from the system. 

Isolated gravitating systems and the structure of their "far fields" have been 
investigated for a long time because of their importance for the interpretation 
of measurements. A series of articles which heavily influenced the way we look 
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at the subject today was published in the early 60's. In these articles vari- 
ous important contributions were made: the "peeling property" of the Weyl 
tensor the idea of analysing the vacuum field equations on outgoing null 
hypersurfaces resulting in the Bondi mass loss formula the invention of 

the NP formalism and the proposal for considering the vacuum Bianchi iden- 
tity as a field equation for the Weyl tensor ||] and the asymptotic solution of 
the Einstein vacuum equations [||. Assuming that certain components of the 
Weyl curvature fall-off in a specific way it was found by formal power series 
analysis of the asymptotic characteristic initial value problem that the fall-off 
behaviour of the fields along null directions could be characterized in terms of 
certain special coordinate systems whose existence was presupposed. Finally, it 
was realized by Penrose [Q that these fall-off conditions as well as the peeling 
property could be understood in a purely geometric way. He introduced the no- 
tion of a conformal extension by which a Lorentz manifold {M, g) is embedded 
into a bigger manifold (M, 17) with isomorphic conformal structure but with a 
Lorentz metric which differs from g by a positive factor g = fl^g. The idea was 
to study the global conformal properties of Minkowski space in order to obtain 
a criterion for what one should call an "asymptotically flat" space-time. Guided 
by the Minkowski situation Penrose suggested that such space-times allow the 
attachment of a conformal boundary J7 which is characterized by the vanishing 
of the conformal factor This boundary is a regular null hypersurface in the 
ambient unphysical manifold. It can be interpreted as the points which are at 
infinity for the physical manifold along null directions. 

The question arises as to what extent this geometric picture is compatible 
with the Einstein equations. Friedrich could derive a system of equations Q, 
the "conformal field equations" , which are defined on that larger unphysical 
manifold. Furthermore, a solution of the conformal field equations gives rise to a 
solution of the standard field equations on the physical space-time. This system 
is written in terms of geometric quantities of the unphysical manifold and the 
conformal factor fl and it is regular everywhere even at points where fl vanishes. 
In a usual 3-1-1 decomposition, the conformal field equations split into constraint 
equations and evolution equations. Using this system of equations Friedrich was 
able to reduce the asymptotic characteristic initial value problem for the Einstein 
equations where data are given on a part of (past) null infinity and an ingoing 
null hypersurface which intersects null infinity in a two-dimensional surface to 
a characteristic initial value problem for a symmetric hyperbolic system Q . 

In order to describe a physical situation one would like to prescribe initial 
data for the conformal field equations on some initial space-like hypersurface 
and determine from them the future of the system. Ideally, the data should 
be given on an asymptotically flat space-like surface. It does turn out that 
the initial data for the conformal field equations on such a hypersurface are 
necessarily singular because the conformal structure of space-time is singular 
at whenever the ADM mass is non-zero (see for a new approach towards 
the solution of this problem). Therefore, the initial data are given on a space- 
like hypersurface which intersects JT" transversely in a two-dimensional surface. 
Such hypersurfaces are called hyperboloidal surfaces because they behave like 
spaces of constant negative curvature in the neighbourhood of J. Friedrich [|lOj 
has shown that the Cauchy problem for data given on such hypersurfaces, the 
hyperboloidal initial value problem, is well posed, i.e., given smooth initial data 
which solve the constraints then there exists a solution of the evolution equations 
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in some neighbourhood of the initial surface. If the data are close enough to 
Minkowski data, then the future development is complete in the sense that there 
exists a regular point whose past light cone coincides with 

The fact that the domain of dependence of the initial hyperboloidal hyper- 
surface includes the complete physical space-time in the future allows the study 
of global phenomena like the behaviour of horizons and the causal structure of 
singularities. But also, since one has access to null infinity where gravitational 
radiation is registered one can in principle "extract" the radiative information 
by purely local manipulations from the fields "on J'". This suggests that the 
hyperboloidal initial value problem is the appropriate device for examining these 
issues. The goal of the present work is the investigation whether the conformal 
field equations, and in particular the hyperboloidal initial value problem, can 
provide an effective numerical tool for analysing the global structure of asymp- 
totically flat space-times and for obtaining information about the gravitational 
radiation emitted by the system in question. The starting point for this investi- 
gation was the work by Hiibner who was able to demonstrate the feasibility 
of that approach in the spherically symmetric case of gravity coupled to a scalar 
field. 

There are various other numerical approaches towards these problems based 
on the numerical solution of the standard field equations, either in terms of a 
Cauchy problem, a characteristic initial value problem or a combination of both. 
The first two alternatives both have some problems. The standard Cauchy 
problem cannot provide complete global information because one has to cut off 
the initial data surface and provide boundary data on a time-like hypersurface 
which intersects the initial surface in a two-dimensional surface. The boundary 
data change the solution in their domain of influence and hence, if the boundary 
data are unphysical, so will be the solution. Even if the boundary conditions 
are physical, the radiation data obtained are still only approximate because the 
boundary is not at infinity. Only there, radiation can unambigously be defined. 
Therefore, as a matter of principle, the standard Cauchy problem can provide 
only approximate radiation information. In the hyperboloidal problem there is 
only one idealization involved namely that of how an isolated system is to be 
described. 

The characteristic initial value problem, on the other hand, can be put to 
good use in the neighbourhood of J. Space-time is foliated by outgoing null 
hypersurfaces and one can perform a conformal transformation to obtain a prob- 
lem where null infinity is at finite places. This is one boundary for the outgoing 
initial null hypersurface. The other boundary is located someplace not too 
far in the interior of the space-time. The numerical procedure for solving the 
characteristic initial value problem is relatively simple compared to the Cauchy 
problem which is due to the fact that the equations split into hypersurface 
equations, which are essentially ordinary differential equations and one evolu- 
tion equation. Three further equations, the so called "conservation equations" , 
have to be satisfied on J'. The problem with this approach is the fact that null 
hypersurfaces invariably tend to form caustics, places where the hypersurface 
intersects itself so that it becomes impossible (or at least very difficult) to give 
unambiguous initial data. The stronger the fields are the earlier the caustics 
will appear. 

The last alternative, commonly called the Cauchy-Characteristic Matching 
(CCM) procedure (probably going back to has been intensely studied 
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by various groups cf. [|T3[ |l^, |T^. The idea is to combine the two previous 
approaches without their respective disadvantages. The procedure is roughly 
to divide the physical space-time by a time-like world tube T and to evolve 
the inner part by solving a Cauchy problem. The exterior of T is evolved by 
solving the characteristic initial value problem based on outgoing null hypersur- 
faces connecting the world tube with null infinity. An initial hypersurface for 
the combined problem consists of a space-like hypersurface iS with a boundary 
(which indicates the intersection of S and T) together with the outgoing null 
hypersurface emanating from the boundary. Obviously, at the interface where 
the initial hypersurface changes its causal character from space-like to null there 
is a non-differentiable kink. The great challenge is to numerically implement 
the information exchange across that kink. 

When viewed in the unphysical space-time the initial surface for the COM 
procedure intersects J" in a two-dimensional "cut". Now consider a space-like 
hypersurface S which also goes through that same cut. It is clear that this 
surface is a hyperboloidal hypersurface. Its domain of dependence is the same 
as that of the Cauchy-characteristic hypersurface. The region of J7 which can 
be described is the same in both cases. One advantage of evolving with the 
conformal field equations is certainly the fact that the causal character of the 
foliation does not change so that there is no interface and no need to change the 
evolution algorithm. Another advantage is that one can go smoothly through 

which allows to keep in the interior of the grid in order to avoid numer- 
ical influences from the grid boundaries. There are more equations to solve in 
the case of the conformal field equations than there are in standard ADM-like 
formulations. This might be considered as a drawback. However, recent for- 
mulations of the Einstein equations as hyperbolic systems introduce many 
additional variables so that the resulting system is comparable in size to the 
system of conformal equations. Furthermore, the quantities in the latter system 
which are evolved in addition to the spatial metric and the extrinsic curvature 
in the standard case have a geometric meaning (Ricci- and Weyl tensor compo- 
nents) and any code which aspires to analyse the space-time structure needs to 
compute those quantities anyway. 

This present article is meant to be the first in a series of papers on the numer- 
ical treatment of the conformal vacuum field equations. In this paper, we derive 
the conformal field equations in a formalism using space spinors. Although this 
has been done previously |l7|| , we present the equations here in a form suitable 
for our immediate purposes, the main reason being to establish a common no- 
tation and for reference. The space spinor formalism has the advantage that 
the equations can easily be decomposed into evolution and constraint parts and 
that the evolution part comes out automatically in symmetric hyperbolic form. 
Furthermore, the equations can be written in a more compact form and the 
possibility of decomposing spinor fields into their irreducible parts can be used 
to remove any redundancy from the set of unknowns. 

In the second article we present the numerical treatment of the evolu- 
tion part with the additional assumption of a symmetry and in the third we 
want to discuss the solution of the constraint equations. The conventions used 
throughout this work are those of Penrose and Rindlcr |l9| . 
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2 The conformal field equations 

In this section we want to give a derivation of the conformal field equations 
and a brief discussion of their properties. Apart from introducing the necessary 
background on the hyperboloidal initial value problem this section also serves 
as a reference to the actual equations used in the code. 

An essential ingredient in this approach towards the examination of global 
properties of space-times is the notion of a conformal transformation. Let (M, g) 
be a Lorentz manifold with vanishing Einstein tensor (we will assume throughout 
that the cosmological constant vanishes). Assume that this "physical space- 
time" is such that the following conditions hold 

• there exists a Lorentz manifold {M, g) with boundary J and a function f2 
on M such that f7 > on M and Vl ^ dVl ^ Q on J , 

• the physical manifold can be identified with the interior of M and there 
the equation g = fl^g holds. 

These conditions state that the physical manifold is conformal to the interior of 
the "unphysical" manifold M. The points on the boundary J' can be thought 
of as representing the points of M which are "at infinity" with respect to the 
physical metric g. With the vanishing of the cosmological constant it follows 
that is a regular null hypersurface in M on which the Weyl curvature vanishes 
(although this is only strictly proven in the case where J has the topology 5*^ x R 



The conformal field equations can now be obtained from the basic geomet- 
ric equations on M and Af , the Einstein equation which holds on M and the 
conformal transformation properties of the geometric fields. In view of the nu- 
merical application it is advantageous to have a first order system. This is easily 
achieved by using a frame formalism. The system then consists of the following 
equations: 

• The first of Cartan's structure equations, which expresses the fact that the 
connection on M is torsion free. It can be viewed as an equation for the 
components of the chosen tetrad with respect to the chosen coordinates. 

• The second structure equation which relates the Ricci rotation coefficients 
of the connection to the curvature components. It can be viewed as an 
equation for the connection components with respect to the chosen tetrad. 

• The Bianchi identity for g. This is an identity which relates the derivatives 
of the physical Ricci and the Weyl curvature. Since M is a vacuum space- 
time this yields an equation for the physical Weyl curvature. Expressing 
this equation in terms of the unphysical connection yields an equation for 
the rescaled Weyl curvature Dated — ^^^Cabcd, which looks formally like 
the familiar spin-2 zero rest-mass equation. 

• The Bianchi identity for g. Again, this is a relation between the derivatives 
of the Ricci and the Weyl curvature, but now on the unphysical space-time. 
Using the equation of the rescaled Weyl curvature this identy yields an 
equation for the unphysical Ricci curvature. 

• equations for the conformal factor Q and its derivatives obtained from the 
conformal transformation law for the Ricci curvature, and 
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• an equation for the function S := jOil^ which is a consequence of the 
earUer equations. 

Due to the geometric origin of these equation, there is gauge freedom in this 
system. Several variables can be chosen freely. Apart from the coordinates, 
this is true also for the tetrad which is fixed by the metric only up to Lorentz 
transformations and for the conformal factor, which is fixed by the conditions 
above up to multiplication with a strictly positive function 0n, where 9 > 
on M. This allows the free choice of eleven functions, the gauge source functions, 
which can be fixed in numerous ways. It is here where the development of a 
code to evolve space-times turns into an art. 

The essential property of the above system is the following: With the gauge 
source functions fixed as arbitrary functions of the coordinates, the system can 
be decomposed by a usual 3 + 1 splitting into two separate systems with re- 
spect to a givem foliation of space-like hypersurfaces. The first of those, the 
constraints, is intrinsic to the space-like hypersurfaces and therefore it restricts 
the values of the variables there. The second part, the evolution equations, can 
be written as a quasi-linear symmetric hyperbolic system. This has the con- 
sequence that the Cauchy problem for this system is well posed: given initial 
data for the unknown functions on a space-like hypersurface S then in a neigh- 
bourhood of S there will exist a solution of the system acquiring the prescribed 
values on <S. It turns out, that once the constraints are satisfied on the initial 
hypersurface they will be satisfied everywhere by virtue of the evolution equa- 
tions, they are propagated by the evolution. Therefore, given initial data which 
satisfy the constraints then they will evolve into a solution of the conformal field 
equations. 

If M is such that the initial hypersurface S and J intersect transversely in 
a regular compact two-dimensional surface then one can talk about the hyper- 
boloidal initial value problem. The standard example for such "hyperboloidal" 
surfaces are the conformal images of the space-like hyperboloids in Minkowski 
space in the usual conformal picture. The fact that S and J intersect in the 
unphysical space-time means in physical terms that the space-like surface ex- 
tends out to null infinity. Thus, such surfaces are not Cauchy surfaces for the 
standard Cauchy problem for the Einstein equations in the physical space-time. 

In summary, the conformal field equations allow a well-posed initial value 
problem on space-like hypersurfaces in the unphysical space-time whose physical 
"pre-images" extend asymptotically towards null infinity. Initial data which 
satisfy the constraints evolve into a solution of the complete system. 

3 Space spinors 

In this section we briefly introduce the basic formalism used to write down the 
equations and to separate them into the evolution and constraint part. This 
can very conveniently be achieved using the space spinor formalism [ pi] , which 
in addition allows writing (and coding) the equations in a more compact form. 

The essential ingredient in the space spinor formalism is a time-like vector 
field which is normalised by tat"' — 2 (note, that we use throughout the 
conventions of |l^). In terms of spinors, we have t" — t^^ and tAA't^^ — 
eA^ ■ The existence of this vector field allows the conversion of all primed 
spinor indices to unprimed ones by extension of the map tta' ^ T^A't^ A to 
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the full spinor algebra. E.g., any covector Va — vaa' is mapped to vab = 
VAA't^ B- Note, that this spinor can be decomposed into irreducible parts: 
vab = ^/"^^abv + vab where vab = vba- In terms of the original covector 
these parts correspond to the components along t"', v = f^Va, and orthogonal to 

to r. 

The vector t^^ can be used to define a complex conjugation map on the 
algebra of imprimed spinors by extension of the map tta ^ t^a '■— tA^ t^A' to 
the full algebra. Note, that ttai...a„ = (— 1)"7''Ai...a„- An even valence spinor 
7rAi...A2„ is called real, iff ttai...A2„ = (-l)"7rAi...A2„- 

The derivative operator V on M can be decomposed as follows 

Vaa' = -^tAA'D — Ia'^ Dab, (3.1) 

or, equivalently 

i^'sV^A' = ^^cabD + Dab, (3.2) 

where D := f^Va and Dab = (bVa)A' are the parts which act along and 
perpendicular to t", respectively. Thus, the general procedure we will follow 
is to write the equations in spinorial form, then convert to space spinors and 
finally decompose them into irreducible parts. 

The derivative of t^"^ gives rise to two important spinor fields, Kab = 
bDIaa' and Kabcd ~ t'^' oDABtcc ■ Note the symmetry and reality prop- 
erties of these fields: Kab = K^^ab) = -Kab and Kabcd = K(^ab){cd) = 
Kabcd- Geometrically, Kab corresponds to the acceleration vector of t" while 
Kabcd is related to the geometry of the distribution defined by ta = 0. This 
distribution is integrable if and only if K^(^bd)a = 0. Then t° is hypersurface 
orthogonal and Kabcd corresponds to the extrinsic curvature of the orthogonal 
surfaces. 

We will assume henceforth that is hypersurface orthogonal. Hence, the 
covector ta is proportional to the conormal of the space-like hypersurfaces given 
by ta = 0. Then, the derivative Dab is the so called Scn-Witten connection 
which plays an important role in various areas of general relativity. It is not 
completely intrinsic to the hypersurfaces but contains information about the 
embedding of the surfaces in space-time. This is reflected in the fact that the 
connection thus defined possesses torsion which is proportional to the extrinsic 
curvature of the surfaces. Therefore, to obtain a completely intrinsic covariant 
derivative operator on the hypersurfaces we define for an arbitrary spinor ttc 
the operator 

QabT^C '■= DabT^C + -j^KABC^T^D- (3.3) 

The connection defined by this derivative operator is torsion free and respects 

the intrinsic metric of the hypersurfaces. Thus, it is the SU{2) spin connection 
of the intrinsic metric. In complete analogy, we define the operator 

d-KC := D-KC + ^Kc^TTD. (3.4) 

This operator defines a connection along the integral curves of which turns 

out to be the spinorial equivalent of the Fermi- Walker connection. E.g., a vector 
ya _ ^AB Fermi- Walker transported along the curves iff dv^^ = 0. 
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These two operators are real in the sense that they map real fields into real 
fields, which is obvious from the relations 

d-Kc = dnc, (3.5) 

dABT^C = "dABTTc- (3.6) 

In order to phrase the structure equations in this formulation we need to 
know about the commutators of these operators because these define the torsion 
and the curvature. The commutators are given by the formula, valid for an 
arbitrary spinor ac 

[d, Oab] ac = ^KabOuc + Kab^^Oefoic - ^abcxc + ^abolc 



QabKcd - QKabcd — Kab Kcdef 



+K^C^KD)EAB + lKABKcD^a'' (3.7) 



and 



2dE(Ad B)01C = ^ABOIC + ^ABOIC 

— ^E(aK^ B)CD — -^KKabCD + -^KefAbK^^ CD 

+ leciAeB)D {KefghK^^''" - K^) | a^, 



(3.8) 



where we have defined the trace of the extrinsic curvature K := K^^ ab and 
introduced the curvature derivations 

^ABac = -^ABC^ctD + 2Kec(A0iB), (3.9) 

^ABac = -^ABC^ao. (3.10) 

Here, we have used the curvature spinors, i.e., the Weyl spinor ^ abcd, the 
(space spinor equivalent of the) Ricci spinor ^abcd '■= ^ABA'B't^ ct^ d and 
the scalar curvature A. 

In view of the reality properties of the derivative operators we find that the 
commutation relation can be split into various parts. This yields equations for 
the spinor fields Kab and Kabcd 

QKabcd = o {QabKcd + QcdKab) + Kab^^ Kefcd + -^^KabKcd 

2 2 (311) 

+ {^ABCD + ^CDAb) - {^ABCD + ABCDj + 4Aec(AeB)D, 

dc{AK^B) = 0, (3.12) 

dE{AK^ B)CD = 2 (^ABCD — ^ABCD^ — ^ i^ABCD — ^CDAb) (3.13) 

and simplified commutation relations 

[d,dAB] OLC = ^KabOuc + KAB^^dEFOic + ]^K:^ {cKd)eabOi'^ 

- ^ {^ABCD - $ABCd) - ^ {^ABCD - ^cdab) (3.14) 
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2dE(Ad^B)ac = 

^-^KKabcd — -^KefabK^^cd — -^^c(a^b)d {K^^'^" Kefgh — ^^)| 

- 2Aec(AeB)Da" + ^ (^abcd + *abcd) + ^ {^abcd + '^cdab) a" 

(3.15) 

When acting on functions, these commutators yield 

[d, dAs] f = \KABdf + KAB^'^dEFf, (3.16) 

^dEiAd'^B)/ = 0. (3.17) 



4 The equations 

We are now in a position to give the derivation of the conformal field equations 
on the unphysical manifold M. To this end we need to introduce coordinates and 
a spin frame with respect to which we express all the spinor fields involved. Since 
we have already assumed the existence of a family of space-like hypersurfaces we 
may now introduce a time coordinate t be requiring that it be constant with non 
vanishing differential on the hypersurfaces. Then, necessarily, we have ta oc Vat. 
Now, we choose arbitrary coordinates {x^,x^,x^} on M. These coordinates can 
be cliaractcriscid by their change along the integral curves of t", thus defining 
lapse function and shift vector. 

The choice of frame is probably best described in terms of orthonormal 
tetrads. On the initial surface wc choose the time-like log of the tetrad to be 
proportional to t"'. Then the other members are tangent to the surface and 
we choose them arbitrarily. To propagate the tetrad off the hypersurface, we 
could use the Formi- Walker transport along the integral curves of f/^ which has 
the property that it leaves the angles between vectors constant, while keeping 
the tangent vector along the curve fixed. However, this is not the most general 
propagation law with these properties. Any other one with the above properties 
differs from the Fermi- Walker transport law by the addition of a term which 
involves an infinitesimal rotation. Thus, we fix an arbitrary transport law which 
leaves angles invariant and fixes the tangent vector to propagate the frame into 
the full space-time. The infinitesimal rotation involved in the transport law 
determines some of the Ricci rotation coefficients of the tetrad thus obtained. 

In terms of spinors this choice of frame is expressed as follows. On the 
initial surface we choose a normalised spinor field oa, i.e., we have o^oa = 1- 
To complete the spin frame we define la = oa- Then the orthonormal tetrad 
constructed from this spin frame has the above properties on the initial surface. 
Now we impose the transport equation in the form 

doA = Fabo"", (4.1) 

where Fab is an arbitrary purely imaginary and symmetric spinor field. It 
corresponds exactly to the infinitesimal rotation mentioned above. It was al- 
ready pointed out that d corresponds to the Fermi-Walker transport, which is, 
consequently, selected by choosing Fab = 0. 
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The connection is defined by the spinor fields Kab and Kabcd, the field 
Fab and finally by the field Tabcd^ defined by 

dABOc = TABCDO^. (4.2) 

It satisfies the relations Tabcd = ^{ab){cd) = —^^abcd, thus being purely 
imaginary. 

The frame components with respect to a coordinate basis are usually ob- 
tained by applying the vector fields which make up the tetrad to the coordi- 
nates. Similarly, in the present case: we apply the operators d and Sab to the 
coordinates 

9t=^, dABt = 0, (4.3) 

dx' = -T\ dABx' = C\b- (4.4) 

This defines several additional fields on M which fix the frame in terms of the 
chosen coordinates. The second equation reflects the fact that we have chosen 
the time-like leg of the tetrad to be the unit normal of the surfaces throughout. 
The function N is the lapse while the three functions T* are closely related to 
the shift vector which appear in all variants of 3 -I- 1 decompositions. 

To make the relationship between the frame as defined above and the coor- 
dinate basis somewhat more precise let us introduce the 1-forms 0, 9^^ which 
are dual to the operators 9, Oab considered as vector fields on M . I.e. for any 
spinor field a"^^ we have the relations 

{6, a^^dAB) = 0, (0-4^, a^^ScD) = a^"" , (4.5) 

{e,d) = l, {e^^,d) = 0. (4.6) 

Then the metric on M when expressed in terms of the present formalism is 
simply 

g = 2e^0-eAB^O'^^. (4.7) 

We have the following relations between the coordinate differentials and the 
forms e, e^^: 

dt = j^e, dx' = C\b^^^ - T"0, (4.8) 

9 = Ndt, 9^^ = D^^ {dx'' + T'^Ndt) , (4.9) 

with D-^^ being the inverse of C^^, i.e., D^^Cab = ^k- I'^om these we get the 
metric expressed in terms of the frame components 

g = 2N^df - Df^DABk [dx' + NT'dt) (g) {dx'' + NT'^dt) (4.10) 

We now come to the system of equations. To express the first structure 
equation in terms of the fields defined above, we apply the commutators to the 
coordinates which gives the equations 

OabN = ^NKab, (4.11) 
diA^KB)c = 0, (4.12) 
d&AB + QabT' = -\kabT' + Kab^^C'ep, (4.13) 



S (aC]5)c = 0. (4.14) 
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The second structure equation is obtained analogously by applying the com- 
mutators to oc yielding the evolution equation for F 

dTABCD — dAsFcD — abe{dF^ c) + -^KabFcd 

+ Kab^^^efcd + -j^K^ (^cKd)eab 

- ^ \^ABCD - ^ABCD^ - ^ {"^ABCD - ^CDAb} (4.15) 

and the constraint equation 

'^dE{A^^ B)CD — '^^^ {B^ C^A)EDF 

+ 2 ^KKabcd — KefabK^^ CD ~ -^^c{a^b)d {K^^'^" Kefgh - 

- 2Aec{A<^B)D + \ {^ABCD + *ASCD } + ^ {'^ABCD + '^CDAb} ■ (4.16) 

The derivation of the equation for the curvature and the conformal factor will 
be given first using primed and unprimed spinors. At the end of this section we 
will collect all the equations in the space spinor formalism. To obtain equations 
for the curvature components we need to look at the Bianchi identities. These 
have the following spinorial form 

Vi3-^*ASCI3 = (b'^CD)A'B- (4.17) 

V^^'$ABA'S' + 3Vaa'A = 0, (4.18) 

which is valid in any space-time. Since the Weyl curvature is conformally in- 
variant the first of these equations when viewed in the physical space-time takes 
the form 

Vb'^*abci3 =0, (4.19) 

with V being the physical connection. Here we have used the Einstein equa- 
tion and the fact that the energy-momentum tensor in the physical space-time 
vanishes. The conformal transformation behaviour of the connection implies 
that this equation is conformally invariant provided we rescale abcd with the 
conformal factor in the appropriate way. Thus, if we define the rescaled Weyl 
spinor 

li'ABCD = n-^-^ABCD (4.20) 

then we have the equation 

^B'^i'ABCD = (4.21) 

holding on A/. Note, that ipABCD is well behaved on M, because the Weyl 
curvature vanishes on J7, where also = 0. 



Expressing the Bianc hi id entities ( 4.17 ) in terms of the rescaled Weyl spinor 



and using the equation (4.21) together with the definition 



'^AA' ■— ^AA'^ 



(4.22) 
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we get an equation for the Ricci spinor 

Vb^'^cda'B' = ^A'^tpEBCD + 2eB(c Vr,)^' A. (4-23) 

When these equations are expressed in terms of space spinors using the operators 
d and Oab they become rather long. To simpHfy them we decompose every 
spinor and every equation into their irreducible parts thus obtaining a set of 
equations which can be further decomposed into real and imaginary parts. The 
resulting equations will be displayed below. 

The next set of equations is the definition of S viewed as an equation for 
Cl and the equation which expresses the conformal transformation behaviour of 
the Ricci spinor viewed as an equation for T,aa' 

VbB'SaA' = -O.^ABA'B' + CabSA'B'S (4.24) 

with S := l/Aon. The final equation is an equation for S which can be derived 
from the equations established so far. It reads 

Vaa'S = -^ABA'B'^''''' + nV^A'A + 2AS^A'. (4.25) 

The conformal transformation behaviour of the scalar curvature implies an al- 
gebraic equation 

2ns -2n^A-EAA'^'^'^' =0- (4.26) 

This completes the derivation of the system of equations and we now want 
to present the full list of equations in the space spinor formalism. The variables 
for which we have an evolution equation are: N, C\g, Kab, Kabcd, ^abcd, 

EaBCD ■= \ [i'ABCD + ''PaBCD^, BabCD ■= ^ {i^ABCD - IpABCD^, (pABCD, 

(pAB, 4>, ^, S, 12ab and S. In this list we have included the irreducible parts of 
^ABCD and T,AA' defined by the decompositions 

SaA' = 2*AA'S — tA'^^AB, (4.27) 
^ABCD = <pABCD + -jf^A{c4>D)B + -^^B{C<pD)A - ■^<^A(C^D)B^- (4.28) 

We have included in this list the lapse N and Kab for which we do not have 
an evolution equation yet. This can easily be achieved by computation of the 
"harmonicity function" F := 20t. Expressing the wave operator in terms of d 
and Oab yields the equation for N, while commuting d and Oab on N gives the 
evolution equation for Kab- 

The evolution equations can be grouped together according to the geometric 
meaning of the variables: 

• The evolution of the frame components 

dN = -KN - N'^F, (4.29) 
dC\B = Kab^'^C'ef - OabT' - ^KabT\ (4.30) 
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• The evolution of the extrinsic curvature and the acceleration vector 

OKab + 2d^'^KABCD = 

Kab^^Kef - KabK - 4<j,AB - 2NdABF - NKabF, (4.31) 



OKabcd - ^ {QabKcd + OcdKab) = 

Kab^^Kefcd + -^KabKcd 
+ 2(f>ABCD - 2Q.EABCD + ^ (6A + (j)) ec{AeB)D, (4.32) 



• The evolution of the intrinsic connection 



1 

2' 

1 

2' 

- iClBABCD - f-{A{c4>D)B) (4.33) 



d^ABCD = QabFcd — '^^ abe(dF^ c) + 7;KabFcd 

+ Kab^^^EFCD + -j^K^ (^cI^D)EAB 

-ifl 

• The evolution of the Ricci curvature 

d(j)ABCD - d(AB4>CD) = 

2 

K(ab4>CD) + K(^AB <I>CD)EF " ■^K(^aBCD)4' 



- Y.EABCD + 2iE(^^BscB)J5, (4.34) 



2 

d(l>AB - -^dABfj) + d'^'^(j)ABCD = 

-Kab4> - 4>abcdK'~'^ + '^Kabcd^'~'^ - -^KcpAB 

+ Eabcd^^'^ - ^dABA, (4.35) 

d(l> + 2dAB<l>^'' = 

- 2Kab(I)^'' - 2K^''^''(j,ABCD + ^Kcj> - 2dA (4.36) 
• The evolution of the Weyl curvature 

dEABCD — 2id(^A^ Bbcd)e = 

2iK(^A^ Bbcd)e — ^K(^AB^^ EcD)EF + "^KEabcd (4.37) 

OBabcd + 2id(^A^ Ebcd)e = 

— 2iK(A^ Ebcd)e — ^K^AB^^ Bcd)ef + "^KBabcd (4.38) 
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• The evolution of the conformal factor 

on = E, (4.39) 

9S = -Kab^^^ + 2S, (4.40) 

d^AB = lKAB^-n<f>AB, (4.41) 

dS =^(t)j:-<l>AB^'^^ +^dA + 2AJ:. (4.42) 



This completes the list of evolution equations. We need to make several remarks. 

• The gauge functions in this system are the "harmonicity" F, the shift 
functions T', the frame rotations Fab and the scalar curvature A. These 
eight functions can be chosen almost at will. From the form of the metric 
we infer a condition which has to be satisfied by the shift functions. Since 
the vector dt needs to be time-like, we find the inequality 

D^^DabiT'T^ > 2. (4.43) 

The full gauge freedom of eleven functions has been reduced to eight be- 
cause of our fixing of the time-like leg of the tetrad. 

• Since the operator d is the directional derivative along the vector field 
t° when acting on functions, most of the equations are simple advection 
equations along that vector field. This is obvious from the explicit form 
of df — jjdt — T'^d^i. There are three subsystems for which this is not 
the case. These are the systems describing the evolution of Kab and 
KabcDi of the Ricci curvature and of the Weyl curvature. They will have 
considerable significance later when it comes to the numerical treatment 
of the equations at the boundary. 

• It is a useful property of the space spinor formalism that the equations 
automatically come out separated into constraints and evolution equations 
and that the evolution equations automatically come out as a symmetric 
hyperbolic system. This is the case for the above system. As written it is 
in symmetric hyperbolic form. The symmetric hyperbolicity is the basic 
property which allows the proof of existence and uniqueness of solutions 
for various initial value problems see e.g. jl^ . 

In a similar way, the constraint equations can be grouped according to their 
geometric meaning starting with 

• the frame components 

= OabN - ]^KabN, (4.44) 
= d(A''C^B)c- (4-45) 

• Next are the extrinsic curvature and the acceleration of the time-like unit 
normal to the surfaces 

= df^A^KB)c, (4.46) 

= d(^A^KB)ECD + I^BabCD - <^{A{C<I^D)B) , (4-47) 
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• the intrinsic connection 

= 2d(^A^TB)ECD + -^KKabCD - -^Kab^^ KcDEF 

'{A \C\^B)EDF 

- ^ec(^es)D {(t> + 6A) + Q.Eabcd + (I^abcd, (4.48) 

• the Ricci curvature components 

= d(^A^(j)BCD)E — 2^(ABC^ 4'D)E — ^BabCD — ^{A^EbCD)E, 

(4.49) 

= O^'^cPaBCD + ldAB(l> + \KABCDcf'' ' ^Kcj^AB 

+ YP'^Eabcd + S^abA, (4.50) 
= a(A^0B)iJ + K(A'^''^cj)B)CDE - i^^'^Babcd. (4.51) 

• and the constraints for the Weyl spinor 

= O^^'Eabcd - iK(A^''''BB)CDE, (4.52) 
= O^'^Babcd + iK^A^''^EB)CDE- (4.53) 

• The rest of the constraints is concerned with the conformal factor and its 
derivatives 

= dAB^ - ^AB, (4.54) 
= dAB^ + KABCD^^'' + fl<l>AB, (4.55) 

= Sab'^cd - ^Kabcd'S + i^<pABCD + ^ec(AeB)D {^<P + 6S) , (4.56) 



= OabS + (pABCD^^^ + ^<pAB^ - ^<t>^AB " ^OabA + 2K^AB- 



(4.57) 



• The final constraint is the algebraic condition mentioned earlier 



= 2OS'-2O2A-is2-SABS^^. (4.58) 

Finally, a note on the name conventions of the various spinor components. In 
the above equations all the spinor fields are irreducible except for the extrinsic 
curvature and the intrinsic three connection which have the decompositions 

Kabcd = KiABCD - -^^c(a^b)dK, (4.59) 

TaBCD = liABCD + 2^^(C72D)B + -j^^B{Cl2D)A — -^^C(A^B)d1 (4.60) 

(4.61) 
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with K4ABCD being totally symmetric. The components of the irreducible parts 
with respect {oa, I'a) are defined as follows for a four index spinor aABCD 

OiABCD = ^4 OaObOcOd " 4q;3 0(^aObOcI^D) 

+ &a2 0(^aobi^ci-d) - 4q;i 0(^a'-bi^ci-d) + «o i-alb'-cI'D, (4.62) 

and for a two index spinor 

aAB = 0:2 oaob - 2ai 0(^a^-b) + "o i^a^b- (4.63) 

When we have spinors with the same kernel symbol but different numbers of 
indices like (f>AB and 4>abcd we specify the number of indices in the name of 
the components, thus obtaining e.g., 4)40 as a component of (jiABCD and 4>22 as 
a component of 4>ab ■ 

5 The symmetry reduction 

Finally we discuss a simplification which has been used to reduce the resource 
requirements. Since the main interests in this project is the development of 
procedures to locate ^ to extract the radiation information from there and 
to study several various gauge choices it is legitimate to assume the existence 
of a continuous symmetry. This is because to locate J as the zero-set of the 
conformal factor is not much more difficult in three dimensions than it is in two. 
On the other hand, one needs to have at least two non trivial spatial dimensions 
because otherwise the Weyl curvature vanishes identically so that there will 
be no radiation present. Hence, in the sequel we will assume the existence 
of a space-like Killing field in the physical space-time, which in addition is 
required to be hypersurface orthogonal. In order to exclude numerical problems 
with coordinate singularities we follow B. G. Schmidt who suggested to look 
at situations where the Killing vector has no singular points. This excludes the 
physically intuitive axisymmetry which has fixed points on the axis and leaves 
us with a completely non physical toy problem. It also has the implication 
that J no longer has spherical cross sections but toroidal ones. These global 
questions are not relevant to local considerations like the influence of different 
choices of gauge functions on the solution or even the question of stability of 
the numerical method. They do, however, forcefully come to the fore when 
it comes to defining and interpreting global quantities like the Bondi mass or 
the radiation flux. These are issues that have not yet been discussed. We will 
consider them in somewhat more detail in a later section. 

We adapt the gauge to the symmetry. First of all, the Killing vector in 
the physical space-time becomes a conformal Killing vector in the unphysical 
space-time. Choosing an appropriate conformal gauge we may achieve that it 
becomes a Killing vector. Then the conformal factor is invariant under the 
symmetry and the conformal gauge freedom is reduced to multiplication with 
functions which are also invariant. We assume that the foliation of M into space- 
like hypersurfaces is compatible with the symmetry, i.e., that is everywhere 
tangent to the hypersurfaces (i.e., orthogonal to t""). Next we choose the frame 
in such a way that one of the space-like legs points along the Killing vector and 
that it is invariant under the action of the symmetry group. Note, that this 
restricts the available frame rotation from the orthogonal group to rotations 
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around the Killing vector. Finally, we take one of the coordinates, say to 
be the coordinate along the integral curves of the Killing vector. Then all 
components of the geometric quantities with respect to the adapted frame are 
independent of that coordinate. 

These choices can be made irrespective of whether the Killing vector is hy- 
persurface orthogonal or not. They do not entail much simplification in the 
equations except for the fact that some frame components and connection coef- 
ficients vanish. In particular, the Weyl curvature still has all of its ten compo- 
nents. However, assuming hypersurface orthogonality simplifies things consid- 
erably. This is because it is equivalent to the existence of a discrete symmetry 
— To explain the simplification it is best to consider an example. Let 
Cabcd be the Weyl tensor on M. The electric part with respect to the Killing 
vector is proportional to ^"'^'^Cabcd which is symmetric under the discrete 
symmetry. However, the magnetic part with respect to which is proportional 
to ^"^'^C'abcd changes sign under the symmetry hence it has to vanish. This 
reduces the Weyl curvature down to five independent components. Similar con- 
sideration can be made for the geometric quantities showing that the number 
of independent variables reduces from fifty-three down to thirty-three. 

In the space spinor formalism we take the Killing vector to be of the form 
= ^^'^ oc o'^^i^\ Then the discrete symmetry implies that the components 
of almost all of the spinor fields vanish if they are obtained by contraction with 
an odd number of o'^ and . Only the fields Fabi ^abcd and Babcd have a 
different behaviour. For them it is the even components which have to vanish. 
So, e.g., the Weyl curvature is described by the five non vanishing components 
Eq, E2, E4, Bi and B3 in agreement with the argument given above. 

6 Conclusion 

The purpose of the present article was the presentation of the conformal field 
equations in the space spinor formalism. This paper is intended to serve as 
a reference for future articles which are intended to discuss the numerical im- 
plementation of this set of equations for solving the hyperboloidal initial value 
problem. We have discussed the simplification obtained from assuming the ex- 
istence of a hypersurface orthogonal Killing vector field. 
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